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1 Introduction 

This paper is devoted to construction of the integral representation for the eigenfunctions of 
the Hamiltonians for the quantum periodic Toda chain. Our work was initiated by the recent 
results of Smirnov [|T| concerning the structure of matrix elements in periodic Toda chain. The 
key point of this approach is Sklyanin's invariant scalar product and we tried to interpret 
it as a part of the Plancherel formula which relates the scalar product of the eigenfunctions 
in original coordinate representation with those arising in the space of the separated variables. 
Such interpretation forces to combine the analytic method of Gutzwiller [^ with the algebraic 
approach of Sklyanin . 

To be more precise, Gutzwiller discovered an interesting phenomenon which understood 
nowadays as quite universal one. If one tries to find the eigenfunctions of the periodic Toda 
chain in terms of a formal series in (auxiliary) eigenfunctions of the open chain with unknown 
coefficients labelled by multiple index, then the multi-dimensional difference equation on the 
coefficients thus obtained can be factorized to the one-dimensional Baxter equations. This 
is what is called the quantum separation of variables. In this way Gutzwiller obtained the 
solution to the eigenvalue problem for the N = 2,3 and 4-particle cases only due to the lack 
of appropriate explicit solution for an arbitrary open Toda chain. The next important step was 
done by Sklyanin [Q . He applied the R-matrix formalism to the present problem and simplified 
drastically the derivation of the Baxter equation to an arbitrary number of particles. The evident 
advantages of Sklyanin's method are a wide area of application and effective way to investigate 
different structures, for example, the invariant scalar products. But, unfortunately, his purely 
algebraic method had no deal with the analytical properties of auxiliary eigenfunctions of the 
specific model and, therefore, some important information is inevitably lost to compare with the 
original Gutzwiller's approach. For example, as a price for efficiency the measure is defined up 
to an arbitrary i^ -periodic function in general and there is no natural way to fix this freedom 
uniquely; there is a fundamental question what is the underlying principle to determine the 
asymptotics and analytical properties of the solution to Baxter equation and so on. Moreover, 
the transition from initial coordinate-dependent wave function to those in separated variables 



is somewhat obscured. The advantage of Gutzwiller approach consists of exphcit information 
concerning the analytical properties of auxiliary eigenfunctions. As a consequence, it gives the 
possibility to calculate many fine tuning quantities, for example the numerical coefficient in 
Baxter equation, the asymptotics of its solutions etc. Therefore, there is a natural problem to 
generalize the Gutzwiller results to an arbitrary A^ -particle Toda chain and to extend them to 
other models. In particular, the investigation of algebraic origin of the auxiliary functions arising 
in other integrable models, is of higher importance. 

We consider the periodic Toda chain as a test example to future generalization in these directions. 
This is the main motivation of the present paper. 

We start from well-known observation that the Whittaker functions 01-0] give the solution of 
the eigenvalue problem for the general open Toda chain [|7|, |^. So it is naturally to use them 
to construct the wave functions for the periodic Toda chain in the spirit of Gutzwiller. As 
the main result we obtain the integral representation for the eigenfunctions of periodic Toda 
chain using intensively the analytical properties of Whittaker functions investigated in [^-0. 
The second important counterpart of this representation is the explicit solution of the Baxter 
equation constructed by Pasquier and Gaudin 0. We calculate the multiple integral by taking 
the residues thus obtaining another representation for the wave function which reproduces the 
particular results of Gutzwiller by putting A^ = 2, 3 and A^ = 4 in our general formula. We 
observe also that the Plancherel formula obtained by Semenov-Tian-Shansky [§] for the wavelets 
of Whittaker functions produce the natural way to calculate Sklyanin's scalar product. It turns 
out that the Sklyanin measure is deeply connected with the Harish- Chandra function arising in 
the theory of Whittaker functions. 

The paper is organized as follows. 

In Section 2 we briefly describe A^ -particle quantum periodic Toda chain in the framework of 
i? -matrix approach and formulate the spectral problem to be solved. The material here is quite 
standard. 

The solution of this problem is described by two theorems stated in Section 3. This are the main 
results of the paper. 

The next four sections are devoted to the building blocks which allows to prove the above 
theorems. 

The solution for the open A^ — 1 -particle Toda chain in terms of Whittaker functions [^ || is 
discussed in Section 4. We introduce Weyl invariant Whittaker functions by simple renormal- 
ization of the original ones and describe their properties using the classical results @]-[|]. We 
calculate also the scalar product of such functions using the Plancherel formula derived in p|. 

In Section 5 we follow the Gutzwiller approach 0] to construct the auxiliary functions which 
play the same role for the periodic Toda chain as the exponentials for the standard Fourier 
transform. We generalize the Gutzwiller solutions to an arbitrary (A^ — 1) -particle solution in 
terms of Weyl invariant Whittaker functions and calculate the action of the diagonal operators of 
the monodromy matrix on the auxiliary functions thus obtaining the functional relations which 
generalize the corresponding relations for the Macdonald function. 

In Section 6 we reformulate the spectral problem using the (generalized) Fourier transform to the 
so-called 7 -representation which is an analytic version of the algebraic approach by Sklyanin 0. 
The analytical properties of the auxiliary function play a crucial role here. We show that the wave 
function in 7 -representation is an entire function with definite asymptotics and it satisfies to 



multi-dimensional Baxter equation. We derive the Plancherel formula thus rigorously obtaining 
the scalar product in 7 -variables introduced in [^. As by-product, we obtain the connection of 
the Harish- Chandra function [^-0 with the Sklyanin measure. 

The solution to the Baxter equation is outlined in Section 7. We essentially use the approach 
by Gutzwiller [^ and especially those by Pasquier and Gaudin 0. 

In Section 8 the proofs of the theorems stated in Section 3 are presented. 

2 Periodic Toda chain (description of the model) 

The quantum A^ -periodic Toda chain is a multi-dimensional eigenvalue problem with A^ mu- 
tually commuting Hamiltonians Hk{xo,po; . . . ; Xn-i,Pn-i) , {k = 1, . . . ,N) where the simplest 
Hamiltonians have the form 

7V-1 

fc=0 
N-l 
H2= Y.PkPrn- J2e->=--^+^ (2-1) 

k<m fc=0 

H3= Yl PkPmPn + --- 
k<m<n 

etc. and the phase variables Xk.Pk satisfy to commutation relations [xk,Pm\ = i^^km ■ This 
system can be nicely described using the i? -matrix approach 0. It is well known that the Lax 
operator 

L„(A) =1^ (2.2) 

satisfies the commutation relations 

R{X - ;u)(L„(A)) ® /)(/ ® L„,(/i)) = (/ ® L„(/i))(L„(A) ® I)R{X - /i) (2.3) 

where 

i?(A) = /®/ + ^P (2.4) 

A 

is a rational i? -matrix. The monodromy matrix 

/ A{X) B{X) \ 
T„,A)^'WA)...MA).(^^^^^^^J (2.5) 

satisfies the analogous equation 

R{X - /i)(T(A) ^I){I^ T{fi)) = (J ® T(/i))(T(A) ® I)R{X - /i) (2.6) 

/^From ( p. 61) it can be easily shown that the trace of the monodromy matrix 

t{X) = A{X) + D{X) (2.7) 



satisfies to commutation relations [t(A),t(/i)] = and is a generating function for the Hamilto- 
nians of the periodic Toda chain: 

N 

tW = E(-l)'^'^"'^^ (2-8) 

fc=0 

Let us consider the spectral problem 

Hk% = Ek% k = l,...,N (2.9) 

where E = [Ei, . . . , E,^) . In other words, \l'g is an eigenfunction of the operator t(A) : 

f(A)vI/^ = t(A;^)vl/^ (2.10) 

where 

t{X;E) = Y.i-l)''X''~'Ek (2.11) 

For the future convenience we shall denote the coordinate dependence of the solutions as 
\I'^(xo,a;) thus selecting xq from all other coordinates x = (xi, . . . ,Xjv-i) • The reason for 
doing this will be clear below. Evidently, the wave function has the following structure: 

'$^{xo,x) = "^^ixo-xi, . . . ,Xn_2-Xm-i) ew\ t^i Y1 ^k \ (2.12) 



h 



k=0 



The main goal is to find the solution to (|2.1CI| ) such that ^£, G L^(R ) . In equivalent terms, 
we impose the requirement 

J /(Ei)^^(xo, x)dE, e L2(M^) (2.13) 

for any smooth function f{y) , (y G M) with finite support. 

3 Main results 

Theorem 3.1 The solution to the spectral problem ( \2.1C\ ), (^.1^ can be represented as the 
integral over real variables 7 = (71, ... , •Jn-i) in the following form: 

^^(xo,x)= ' _\., J fi-\l)C{r,E)^^,EA^o,x)d'y (3.1) 

where 
(i) The function /i(7) is defined by the formula 



Ki) = n 

j<k 



ih 



(3.2) 



(ii) The function C^j^E) is the solution of multi- dimensional Baxter equation in the Pas quier- 
Gaudin form J^ 

^^' ' ^ ^^ ~^^. Vf r^^ (3-3) 



i=i n^sinh|(7,-5fc(£;; 



where the entire functions c±{'')) are two Gutzwiller solutions ^ of one- dimensional Baxter 
equation 

t{r, E)c{r, E) = r^c(7 + ih; e) + i^cij - in-, e) (3.4) 

and the parameters C,{E), 6 = {6i{E), . . . ,6,^{E) ) satisfy the Gutzwiller conditions (the 
energy quantization) /^, ^. 

(Hi) The wave function \l/-y^£;^(a;o, a;) is defined in terms of Whittaker function w{x;^) /0/-/|^ 
for the GL(A^ — 1;R) open Toda chain /0, ||/ 

_ 1 • N-l 

j<k mi n ^^^ 

The function (^.d^ ) is symmetric under the permutation of 7 -variables Q. 



The multiple integral ( |3Tl| ) can be explicitely evaluated. Let y = {yi, . . . ,yi^) be an arbitrary 
vector in R^ . We denote y^'') = (yi, . . . , y^_i, y^+i, ...,y^)e R^'^ 



Theorem 3.2 Assuming that Sj{E) 7^ Sk{E) , the solution ^3. 1\ ) can he written (up to unessen- 
tial numerical factor) in equivalent form 

%{xo.x) = Y.{-^r-' E A(5W+^;.nW)C+((5W+zW^))vl/,(.,^^^^(,,_^^(xo,a.) (3.6) 
s=i Ti(=)ez^-i 

where 

N-l 

C4j)^ llc+iir.E) (3.7) 

i=i 

and A (7) = JJ ['-fj — 7^) is the Vandermonde determinant. 
j>k 



Remark 3.1 For N = 2,3 and N = A particle Toda chain the formula ( ^.(Jj ) reproduces the 
results obtained by Gutzwiller J^. 



N-l 

1 In (|U) we use the notation (7,p) = 5 E {N~2k)-ik ■ 

fc=i 



4 ^^(A^ — 1,m) Toda chain and Whittaker function 

It can be easily shown that the operator 

Af-l 

5(A) = e-^o ^(-l)^A^-'=-i/i,(a;i,pi;...;a;^_„p^_J (4.1) 

fe=0 

is the generating function for the Hamiltonians hk of GL{N — 1) Toda chain 0. In particular, 
the simplest Hamiltonians have the form 

N-l 

4.2 

N-2 ^ ' 

j<k fc=l 

Let 7 = (71, . . . ,7iv-i) be the set of real parameters. We consider the spectral problem for 
GL{N -l,R) Toda chain 

hkiJ~,{x) = akh)i^^ix) k = l,...,N-l (4.3) 

where crkil) are elementary symmetric functions. In equivalent terms, '?/'t(^) satisfies to the 
equation 

Af-l 

BiX)i;^{x) = e-"° Y[i\- 7^) ij^{x) (4.4) 

m=l 

As a trivial corollary, 

5(7,)V^^(x) = V7, E 7 (4.5) 



Remark 4.1 The equations ( U-^j ) are an analog of the notion of the "operator zeros" introduced 
by Sklyanin j^. 



Obviously, in the asymptotic region x^+i ^ x^, {k = 1, . . . , N—2) all potentials vanish and the 
solution to (^4.3|) is a superposition of plane waves. The main goal is to find the solution to ( [4.3|) 
satisfying the following properties: 

(i) The solution vanishes very rapidly as Xk — Xk+i — > 00 for any given k . 
(ii) The function ip-y is Weyl-invariant, i.e. it is symmetric under any permutation 

V'...7,...7fc... = ^-^k-nj- (4-6) 
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The commutativity of the operators hk follows from the commutation relation [B{X), B{fi)] ~ which is 



encoded in (2.6) 



•^ More precisely, ^~^{x) ^ cxp{—^e'-^'' ^'»+i)/2} as x^ ~ x^+i 



The property (i) defines the unique solution (up to common 7 -dependent factor). It is called 
the Whittaker function w{x] 7) whose analytical and invariant properties were studied in [^-0- 
In the context of the open Toda chain the Whittaker functions have been appeared for the first 
time in 0, ^ . The requirement of Weyl invariance ( [4.6|) is very convenient for our purposes. It 



turns out that symmetric function ip~^ is most suitable to construct the solution for the periodic 
Toda chain. 

For GL{N —l^M.) Toda chain the Whittaker function has the following integral representation 

i-iB 

7V-2 



w{x-n) = e^"^ / rfz n Af^^-^'=-)-^(z) exp{|;.,,,+ie(^-^-^)/^} (4.7) 



where the integration is taken over the upper triangular (A^ — 1) x (A^ — 1) matrix z = \\zjk\\ 
with unit diagonal and Afc(z) are the principal minors of the matrix zz^ . Finally, (., .) is the 
standard scalar product in M ~^ . 

In particular, taking A^ = 3 and using the integral representation for the Macdonald function 
Kyiij) , one obtains the appropriate solution 

71 -72 

w{x- 7) = 2v^ , ^ ^^, e^(^^+^^)("^+"^)ir7i-72 (le^-^-^y^) (4.8) 

p/71772 _^ i\ -^n-^h ' ^ ' 

of the spectral problem 

(Pl + V2)W = (71 + 12)W 

r _ 1 (4-9) 

{P1P2 -e""' ^'|w = 7i72W 

Indeed, Ky{y) ~ y^^/'^e'^ as y ^ 00 and property (i) is fulfilled. 

All essential theorems for the classical Whittaker function 0]-[^ hold in the present case as well. 
We represent four lemmas below without proof 0. 

Lemma 4.1 For any s &W 

w{x;sj)=M{s;j)w{x;j) (4-10) 

where for the permutation Sjk '■ 7j <-> 7fc (j < k) 



k r ( Tj ~T"' _|_ lArf 7m-7fc I i' 

M{s,,;-f)= n n'^^^"^-^^^^ ) '' ^2 '^ ^^J (4.11) 

m=3+l '- [~^f;~ + 2 ) '^ \~^:h~ + 2 , 



^ The function (4.7) differs from the corresponding solution for the SL{N —IjM) case by the factor 

JV-l JV-l 






The notations and results are extracted from the papers 



The coefficients A4{s]j) satisfy to relations 

M{siS2] 7) = M{si; S2-f)M{s2] 7) (4.12) 

M{s; -f)M{s; 7) = 1 Im7fe = (4.13) 

(where z always denote the complex conjugation of z ). ■ 



Lemma 4.2 The Whittaker function w{x; 7) can be analytically continued to an entire function 
of IE C"-' . ■ 



Lemma 4.3 There is a unique solution v{x] 7) to ( U-^j ) with the asymptotics 

in the region x^+i ^ Xk, {k = 1, . . . ,N — 2) . Let W be the Weyl group (the permutation 
group of the variables 'y ). The functions v{x; 37), (s G W) form a basis of the solutions to the 



spectral problem (4-5) 



Lemma 4.4 The Whittaker function ( H-lD has the following expansion in terms of basis func- 
tions v{x] sj) , (s G W) : 



w{x--i) = Y. c{s-i)M ^(s;7)t;(a;;s7) ui^\ 

seVF 

where 



cW^n^^.i) (4.16) 



is (non-normalized) Harish- Chandra function^. 



Now we are able to define the Weyl symmetric solution and to describe its properties. Let p is 



a half-sum of positive roots of sl{N — l,M.) written in standard basis {e^} G 



oN-l 



p=J2p,ek; pk = -{N-2k) (4.17) 

fc=i ^ 



Lemma 4.5 The function 

^^{x) = ;,-2^(^.p)/^ n ^"'^'r(^+i) w{x-n) (4.18) 

j<k 
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In (|4.16| ) B{z, 



,w 



- r(z)r(^;) 
~ r{z + w) 



(i) is Weyl invariant. 

(a) it can be analytically continued to an entire function of j E C ~^ . 
(Hi) in terms of the basis v{x; 37) , {s G W) 

V^^(^) = H ^(S7)«(s7)^(»;;s7) (4.19) 

sew 

where N{j) = h-^'^'^'^^^^ and 

«(7) = nr(^) (4.20) 

j<k 

(iv) for real 7 variables the function ip-^ obeys the asymptotics 

2 — N (IT '\ 

^^{x) ~ |7,|— exp{ - _(iV-2)|7fc|} (4.21) 

as |7fc| — »• 00 . 



Proof. The Weyl invariance is a direct corollary of the formulae ( |4.10| ), ( [4.11| ). 



To prove (ii), one should note that the poles coming from F functions in the r.h.s. of ( [4.18D are 
cancelled by the corresponding zeros of the Whittaker function w{x] 7) . Indeed, the function 
A4(sjfc;7) has the poles at 7^ — 7a: = —ih{n + ^) , {n E N) . But the function w{x]Sjk'y) is 
an entire function according to Lemma ^^. Hence, from (|4. lOj) it is clear that w{x] 7) vanishes 
exactly at the same points. 
To prove (iii), let 

6(7)-n--^/^r(^ + l) (4.22) 

j<k 

for brevity. Then the Harish- Chandra function ([4.1 6|) is written identically as a ratio 



while the definition ( [4. 18] ) reads 

^Ij^{x) = N{j)b{j)w{x; 7) (4.24) 



Using the Weyl invariance ips-y = i^-y , one easily finds from ( [4.10 ) and 



-'-)^i;^i^ 



Hence, the expansion ( [4.19[ ) is a simple corollary of ( [4.15[ ) and definition ( [4.24[ ). 



The proof of (iv) follows from the Stirling formula for F -functions and from the calculation of 
the integrals in ( [1.7[ ) by the stationary phase method []. ■ 



^Actually, one can prove the exponential decreasing ^-y ~ exp <— 4t(A^ — 2)|Rc7fe| > as |Re7fe| ^ 00 in any 
finite strip of the complex plane. 



Lemma 4.6 Let 7 = (71, . . . , 7jv-i) and 7' = (7J, . . . , 7^_ J be the real parameters. Then 

J ^y{x)^^{x)dx = {2nhf~^\a{-f)\^ ^ 6{s-f - 7') (4.26) 



sew 



Proof. This is a consequence of the Plancherel formula proved in p[ for the SL{n,M.) case. 
The scalar product (^4.26|) results from Semenov-Tian-Shansky formula taking into account the 
relations ( ^T2|) , ( p^ and (^^). ■ 



Conjecture 4.1 The Whittaker functions obey the completeness condition^ 

00 
J \a{j)\-^ij^{x)i^^{y)d-f = (iV-1)! {2nhf-'6{x - y) (4.27) 

—00 

5 Auxiliary functions 

Let us introduce the auxiliary function 

JV-l 
T / N dcf i^°(^-E7n.) (5.1) 

where e is an arbitrary real parameter. Obviously, this function is the solutions of ( [4.3|) or, in 
equivalent terms, 

N-l 

5(A)^^,,(xo, x) = e-"° n (^ - 7m) ^^..(a^o, a;) (5.2) 

m=l 

In particular, 

E(7,)^-,,,(a;o,a;) = V7, G 7 (5-3) 

Further, by construction the function ^-^^t satisfies to the equation 

iJi^^,, = e^^,, (5.4) 

The scalar product of the auxihary functions follows from (f4.26 ): 



j ^Y,e'(a;o, aj)^-^,.(xo, x)dxodx = {2nh)^ n{j) 6{e-e) J2 ^i^l - V) (5.5) 



_ sew 

where 



■p / 7i 7a; ' 



Moreover, the auxiliary functions satisfy to completeness condition 



(5.6) 



I fi-\j)^^,,ixo,x)%,,{yo,y) dyde = (N-l)] {2TThf5{xo - yo)5ix - y) (5.7) 



® We confirm this conjecture by explicit calculations for GL{2, M) chain. The general proof will be published 
elsewhere. 
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provided the conjecture ( t4.27[ ) holds. 

Certainly, the function (5J.) does not satisfy the whole system ( 2.10 ). Nevertheless, it possesses 



some remarkable properties which allow to construct the eigenf unctions of the periodic Toda 
chain. 

First of all, the functions "^-y^e and "^-y^e can be extended to the entire functions of the complex 
variables 7 G C ^^ . This is an evident consequence of Lemma [4.5|. We shall denote the 



corresponding analytic continuations by the same letters. Further, the following statement will 
be of importance below: 



Lemma 5.1 The action of the diagonal operators of the monodromy matrix l \2.l!^ on auxiliary 
function ^-^^e (7 £ C^^"*^) is defined by the formulae 

A(7j)^^,, = i^^^+ine^^e (5.8a) 

D(7,)^^,. = r^^^_.,;,,^,, (5.8b) 

Proof. We consider the following commutation relations encoded in ( p.6| ): 

(A - /i + ih)D{ij)B{X) = (A - iJ.)B{X)D{fi) + inD{X)B{id) (5.9) 

(A - /i + ih)A{X)B{ij,) = (A - i2)B{ij)A{X) + iM(/i)5(A) (5.10) 

Taking /z = 7^ in (|5.9| ) and applying this operator identity to the function ip-y{x) , one arrives 
to the relation 

B{X)Di^,)i;^ix) = e-^« (A - 7, + th) J] (A - Im) D{^^)ij^{x) ^^_^^^ 

since ( |4.5|) holds . Hence, D{jj)ip~^{x) is a solution to the spectral problem (^4.4|) with the set 
of eigenvalues 

7 -^ (7i,...,7j_i,7j-i/i,7j+i,...,7jv-i) (5.12) 

i.e. one obtains the expansion 

D{,lj)^~i{,x) = Y, dj{s; 7) v{x; 5(7 - ihej)) /g^^gx 

sew 

due to completeness of the basis v{x; 7) (s G W) . On the other hand, the function D{'^j)ip~^{x) 
vanishes as Xk—Xk+i — ^ C)0 for any given k . Hence, the r.h.s. in (|5.13|) is proportional to ip~^^ijiej 
because of uniqueness of the Whittaker function. Thus, 

D{-^j)il)^{x) = djij) ip-^-ine, (5.14) 

To find the coefficients dj{^) , we consider the asymptotics in the region Xk+i ^ x^ , {k = 
1, ... A^ - 2) where ( |4.14| ) holds. Note that in this limit 



'N-2 



D{X) = -e"'^-!-"" <^ n (A -P-) + C'(max{e"^-"'=+i}^^J (5.15) 



^m=l 



^Note that D{^) does not contain the operator po and, therefore, commutes with e ^° . 
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Let 7 be the complex variables with Ini7A: < . Then it is easy to find that asymptotically 

X''A'7l,...,7j_l7j+l,...,7jv-il^l! ■ ■ ■ ,Xj, 



mj^j ~ ^^ ' (5.16) 



where ijj is the asymptotical value of GL{N — 2,W) Whittaker function. On the other hand, 
under above conditions 

mi ^ '^ ^ (5.17) 

X'?/-'7i,...,7j_i7j_,_i,...,7jv^i(,2J1; • • • ,Xff_2) 

The comparison of ( p.l6| ) and ( |5.17| ) gives ^^(7) = i'^e'^" . Hence, in the region Im7fc < 

D{^,)tP^ix) = t-^'e-'^y^.in^^ix) (5.18) 



and this relation can be obviously continued to the whole complex region. The formula (|5.8b|) 
follows then from the definition ( |5.1| ). 



Of course, it is possible to prove (|5.8a|) according to the same reasoning. Actually, the relation 



(|5.8a|) is a corollary of (|5.8b|) since the quantum determinant of the monodromy matrix (|2.5|) is 



equal to 1. Applying the operator identity 

A(7, - th)Di^,) - C(7, - ih)Bi^,) = 1 (5.19) 

to the function \E'^ ,= and using the relations ( |5.3| ) and ( p.8b| ), one arrives to ( |5.8a| ). ■ 



As a direct consequence of ( |5.8|) and definition ( p.7|) , for any 7 G C 



^N-l 



' ' (5.20) 



f^(7j)^^,, = i^<i/^+ine„e + Z"^^^- 



where t^{X;xo,po;...;Xj„_i,pN-i) = f(A; Xq, — Po; • • • ; a^iv-i, — Piv-i) 5 (-^ ^ C) and the analytic 
continuation of the functions is understood. 

Remark 5.1 The formula ^5.2C\) has been derived by Gutzwiller ^ using cumbersome explicit 
calculations for n = 2 and n = 3 -particle open Toda chain. 



To prove Theorem ^]T| one needs to determine the action of the operator t(A) (A G C) on the 
function "^-y^e ■ 

Lemma 5.2 

Af-1 N-1 N-1 \_ (5 21) 

= (A-e+E7.) n(A-7™)^-,,e+E n ^^0^^^+..e,, + r^Vl>^-..e, ' 
j=l m=l j=l m^j il '"^ 
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Proof. The function t(A)^-ye is a polynomial in A of order N with the leading terms 

f(A)^^,, = (A^ - A^-^e)^^,, + 0(A^-2) (5.22) 

according to (^.41). Moreover, this polynomial is defined at the points A = 7^ , (j = 1, . . . , A^ — 1) 
due to ( ^.201) . Using the standard interpolation formula, one easily obtains (|5.21| ). ■ 



Remark 5.2 The formula ( \5. 2]\ ) is a coordinate version of the operator interpolation formula 
obtained in Ml- 



6 Periodic Toda chain in 7 -representation 



Let \l/g, {xq^x) be the fast decreasing solution of the problem ( p.lO|) . We define the function 
C(7) by generalized Fourier transform: 

00 
5{Ei - e) C{-f] E)= j *g,(xo, x) ^^,e(xo, x)dxodx (6.1) 



Theorem 6.1 The function C{'y) possesses the following properties: 

(i) It is a symmetric function with respect to 7 -variables, 
(a) It is an entire function of the variables 7 G C^"^ . 
(Hi) The function C(7) obeys the asymptotics 

Ch) ~ |7fcr^/V^^ (6.2) 

OS |7fc| — > 00 along the real axes, 
(iv) The function C{^) satisfies the multi- dimensional Baxter equation 

t{^j- E)C{r, E) = i^Ci-j + ihej- E) + r^C(7 - iUej- E) (6.3) 

where ^(7; E) is defined by ( \2.1]\ ). 



Proof. The symmetricity of the function C(7) is obvious. 

We present here only the sketch of the proof concerning the statements (ii) and (iii). 
The statement (ii) follows from the assertion that the auxiliary function ^-y^e is an entire one 
while the solution of the periodic chain vanishes very rapidly as \xk — Xk+i\ — > 00 . [^ 
(iii) The asymptotics (|6.2| ) is a combination of two factors. The first one comes from the asymp- 
totics ( [4.21|) which results from the Stirling formula for the F -functions in the definition ( [4.18| ). 



^° Actually, the boundary conditions have the same importance here as the requirement of compact support 
in the theory of the analytic continuation for the usual Fourier transform. 
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The second factor ~ |7fc| ^ exp{— 7r|7fc|//i} results from the stationary phase method while cal- 



culating the multiple integral including the Whittaker function ( |4.7| ). The calculation is based 
heavily on the exact asymptotics of the function \E'g (xq, x) as \xk — Xk+i\ — > c>o . 
The proof of (iv) is simple. Using the definition ( p.lO|) and integrating by parts (evidently, 
boundary terms vanish), one obtains 



(6.4) 



6{Ei - e) t(7j; E)C{-f) = J {t{-fj)^^{xo, x)} ^^,e(xo, x)dxQdx 

— oo 

*g(xo, x) {r(7j)^^^,(xo, x)}dxodx 



Taking into account the relation (|5.2CI|) , the Baxter equation ( |6.3|) follows from definition (|0 



Let us consider the following function [^ 



%{xo,x) = j^^ J f^ \-f)CiT,E)^~,,EAxo,x)d-f 



The integral ( |6.5| ) integral is correctly defined. Indeed, the measure 



(6,5) 



.-'(7) = nr(^) 



j<k 



ih 



TT 



(6.6) 



(7r/i)-(--)<--'/^n (7.-7.) sinh ^(7,-7.) 
j<k " 



is an entire function. Therefore, there are no poles in the integrand. Moreover, 

/^-'(7)~l7.r-'exp{|(iV-2)|7.|} 



(6.7) 



as |7fc| ^ oo . Taking into account the asymptotics ( [4.21| ) and (|6.2| ) one concludes that the 
integrand has the behavior ~ |7fc|^-'^exp{— 7r|7fc|/^} as |7fc| — > oo . Therefore, the integral ( |6.5D 
is convergent. 

Using the scalar product ( |5.5| ), one can write the Plancherel formula Q 

oo oo 

{2nhf J %,{xo, x)'^Jxo, x)dxodx = — — 6{E,-E[) J ^-^(7) C{t, E')C{t, E)d-f (6.8) 

—00 —00 

The integral in the r.h.s. of ( |6.8| ) is absolutely convergent due to asymptotics ( |6.2[ ) and (|6.7|). 
Hence, the norm H^E'^.H is finite modulo GL{1) 5 -function 6{Ei—E[) (see corresponding factor 
in ( p.l2|) which results to this function). 



^^The expression ( |6.5| ) is an inversion formula t o (p.ll) assuming the completeness condition ( p.7| ). Otherwise 
one can only conclude that (6.1) is a corollary of (6.5) due to the scalar product (^.5[) 



^^ Let us stress again that (6.8) is deeply connected with the Plancherel formula for SL{n,W) Toda chain 
derived in IHl. 



14 



Remark 6.1 The scalar product in 'y -variables appears for the first time in J^. We see that 
the Sklyanin measure fJ^^^i"/) is naturally connected with the Harish- Chandra function ^4-2^ . 



To prove that the function (|6.5|) satisfies to the spectral problem ( p.lO|) , one needs more details 
information concerning the asymptotics of function C(7; E) : it should exponentially decrease 
in any finite strip of the complex plane 7^ G C as | Re 7^ | — »• 00 . Then ( |6.5|) is a solution if 
C(7; E) satisfies the Baxter equation ( |6.3| ) (see the proof of Theorem |3.1| below). Up to now 
we can prove only the asymptotics (|6.2|) starting from the definition ( |6.1| ) (though the analytic 
continuation seems almost evident). Therefore, in the next section we describe the explicit 
solution to equation ( |6.3| ) p, |^ which obviously possesses the required asymptotics. 



7 Solution of Baxter equation 



In this section the Pasquier-Gaudin solution to the Baxter equation ( |6.3| ) is described which 
leads to refined derivation of Gutzwiller quantization condition 0] . 



We solve the equation (|6.3|) in the separated form 



N~l 

CiT,E)= nc(7.S^) 



(7.1) 



Obviously, the "one particle" functions 0(7) satisfy to the one-dimensional Baxter equation 

t(7; E) c(7; E) = z^c(7 + zh; E) + r^c(7 - ih- E) (7.2) 



N 



where ^(7; E) is, by definition, a polynomial H (7 — ^k{E)) . We assume that the roots Xk{E) 

k=l 

are real for the solutions with the finite norm 0. We shall oftenly do not write the explicit 
dependence on energies for brevity. 



Lemma 7.1 /j^ The equation l \7.2{ ) admits two fundamental solutions 



c±(7) 



ttN^ 



c±(7) 



where 



c±(7) 



K±h) 



N ^n. 



llK'^T{lTJ^il-\k] 



k=l 



(7.3) 



(7.4) 



and K±{'~f) be the following N x N determinants: 

1 



^+(7) 



1 

1 


t{-f+ih) 

1 
1 




1 


n 


t{-f+2ih) 



t{^+2ih) 
1 


1 



t{-f + 3ih) 



t{-f + 3ih) 



(7.5a) 
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This can be demonstrated directly for A^ = 2, 3 particle chain. 
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^-(7) = 



t(7- 


-3ih) 





1 
1 


t(7 


1 




U 

1 






til 


-2th) 


t(7 


-2 


i/i) 













1 




1 





t{-i-ih) 



(7.5b) 



Proof. The determinants -^'±(7) are correctly defined |jTO[ and satisfy to the equations 

1 



K^il-th) = K^ii)- 
K_{^ + th) = K_{^)- 



t{^)t{^+th) 

1 
t{^)t{-f-tn) 



K+{-i + ih) 
K^{-f-ih) 



(7.6) 



This can be easily proved by expansion of these determinant with respect to the first row. The 
rest is trivial. Note that solutions ( [7.4| ) are determined up to arbitrary z^ -periodic function. ■ 



(7.7) 



Remark 7.1 The function c±{'~f) are the solutions of the Baxter equation 

t(7) c(7) = r^c(7 + ih) + i^c{-f - ih) 

Lemma 7.2 The functions 0^(7) possess the following properties: 

(i) C-|-(7) are the entire functions in 7. 
(a) Let 7 he real. Then 



c±{l) 



7rJV|7| 

-N/2 ^^T- exp 



7rAr|7| 

|^|-A^/2g 2n exp 



iN'y I7I niN 
iN'-f I7I niN 



7 



7 



+ CX) 



—00 



(7.^ 



(Hi) 



N\l\ 



^log|7l| 

c_(7) ~ exp| — log|7|j 



7 -^ +100 
7 — > —ioo 



(7.9) 



Proof, (i) It is easy to see that the poles of -^'±(7) are cancelled by the corresponding poles 

coming from the product of F -functions in (|7.4|) . Hence, C-|-(7) are an entire functions. 

(ii) It is clear that lim -^'±(7) = 1 for real 7. Similarly, -^'±(7) ^1 as 7 -^ ±zoo . The 

7|— >oo 

asymptotics follow from the Stirling formula for F -functions. Note that the Stirling formula gives 
even more detailed asymptotics than (|7.8|) ; actually c±(7) ~ exp{^|Re7|} as |Re7| ^ 00 in 
any finite strip |Im7| < e, (e > 0) . ■ 
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In order to construct the solution with the asymptotics ( |6.2D , one should recall that the solutions 
([77^ ) are defined up to any z/i -periodic factor. Clearly, the factor 



7r7 



N 



(7.10) 



(where 5k are an arbitrary parameters) does not spoil the asymptotics of the functions (|7.4| ) as 

27r|7 

7 -^ +00 along the real axis while giving the correcting factor e ^ as 7 ^ —00 . Hence, 
one can consider the following solution of the Baxter equation (|7.2| ): 

c+(7) -^c_(7) 

(7.11) 



c(7) 



nsinh|(7-5, 
fc=i "' 



where ^ is a arbitrary constant. Due to ( [7.8| ) this solution has prescribed asymptotics 

7rA''|7 

c(7) ~ |^|-^/2e-^^ 



(7.12) 



as 7 -^ ±00 . On the other hand, the denominator in ( [7.11|) has the infinite number of poles at 
1 = 5k + iTiUk , nfc G Z, /c = 1, . . . , A^ and the solution ( [7.11[ ) is not an entire function in general. 
The poles are cancelled only if the following conditions hold: 



c+{5k + ihuk) = ic^{5k + ihuk) 
In turn, this means that the Wronskian 

1^(7) = c+(7)c_(7 + ih) - c+(7 + ih)c^{'y) 
vanishes at 'j = Sk + ihuk ■ Due to Baxter equation ( [7.7| ) 

w{-f + ih) = (-1)^1^(7) 

Substitution of O to (|7.14|) gives 



TV 



TT 



W{-f) = i-^Uin) TT T^'^h sinh -(7 - \k 

k=l ^ 



where the function 



n{j) = K+{-f)K_{j + ih) 



K+{^ + ih)K_{^) 



is called the Hill determinant ||1CI|| . 

Lemma 7.3 /|7^/,/|^/ The function 7i{'j) has the following three- diagonal determinant f( 



(7.13) 
(7.14) 
(7.15) 

(7.16) 
(7.17) 



orm: 



nil) 



1 


1 






t{^-ih) 




1 


^ til) 





n 


, V, 1 


1 



ti-f+ih) 



t{-f+ih) 



(7.18) 
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This infinite determinant is absolutely convergent; 'H{'-f) is a meromorphic function with the 
poles at the points X^ + in^h , {uk G Z, A; = 1, . . . , A^) . The Hill determinant is ih -periodic 
functions with the asymptotics 'H{'^) = 1 + 0(7""'^) as 7 — i> ±00 . The following representation 
is valid: 

N 

n{^) = l + J2ak{E) coth-(7-Afc(^)) (7.19) 

fe=i 

where ak{E) are some definite functions depending on the parameters E only. These functions 
satisfy to the constraint 

N 

E«fc(^) = (7-20) 

fc=i 

The equation 7^(7) = has exactly N real solutions Sk{E), {k = 1,...,N) subjected to the 
constraint 

N N 

J25,{E) = J2X,{E) (7.21) 

fc=i fe=i 



Due to its periodicity, 7-^(7) = at any point Sk{E) +ihnk {n^ G Z) . Hence, the solution ( [7.11[ ) 
has no poles if one takes 6^ = Sk{E) provided that the constant ^ is chosen in such a way that 



t _ c+(^) 



c_(7) 

Hence, one arrives to the following 
Lemma 7.4 J^ The function 



k = l,...,N (7.22) 



^, p. fr c^{lfE)-aE)c.{^,-E) 

^^^' ^> = 11 -^ TT — (7.23) 

^=1 nsinh|(7,-4CE^ 



where 5k{E) are the zeros of the Hill determinant ([7. i<^ and the constant ^ is chosen according 



to (\7.22J , satisfies to conditions of Theorem ^A 



The quantization conditions 

c_(5i) ■■■ c_(5,) ^ • ' 

determine the energy spectrum of the problem. These have been obtained for the first time by 
Gutzwiller H using quite different method. 



8 Proof of the main theorems 



Proof of Theorem |3.1| . We have proved already that the function "^^ is integrable in a sense 
of ( |6.(j| ), ( p^.l3| ). Now we prove that the function ( |6.5| ) satisfies to the spectral problem ( |2.10| ) if the 
entire solution C{-j; E) to the Baxter equation ( |6.3|) obeys the asymptotics ~ exp{ — ^|Re7j|} 



as |Re7j| — > cxD in the strip —iTi < Im7j < ih . Note that the Pasquier-Gaudin solution ( [7. 231) 
satisfies to this requirement (see proof of Lemma [77 



Due to (|5.21|) one obtains 

CX) 

. -1) 

Af-l N-1 N-l X—'y 

j=l m=l j=l m^j ij '^"^ 

Changing 7^ -^ 'jj ± ih in appropriate parts of the integrand and noting that the measure 
satisfies to the difference equation 

^.-\^ + ^n5,) = (-l)V^(7) n \~!"1^'^ (8.2) 

one arrives to the expression 

°? Af-l Af-l 

f(A)^^= y rf7;,-i(7)C(7)(A-Ei+E7,) n(^-7j^7,i^.+ 



-CX) 



j=l m=l 



+z^ 5: / rf7i • • • / ^7i • • • / ^7iv-i /i-^(7)C(7+^^e,)^^,s, H ^^+ (8.3) 

1=1 ,- 1713^1 7i 7m 

iY_i 00 oo+i?i 00 

+z-^ E / ^71 • • • / rf7, . • • / ^7.-1 fi-\-f)C{-f-the,)^^^E, n ^^ 

7=1 , -f m^i 7j 7m 

The integration contours [—00 ± i/i, 00 ± i^] can be deformed to the standard one [—00, 00] . 
Indeed, the corresponding integrands are the entire functions (all possible poles are cancelled by 
appropriate zeros of the measure /i~^(7) , see ( |6.6|) ). Moreover, it is not hard to see that the 
integrand is exponentially decreasing in the strip —ifi < Im7j <ih as |Re7j| — ;> 00 . Therefore, 
the contour integral around the strip vanishes and all integrals in (|8.3|) can be deformed to the 
standard ones. Finally, using the Baxter equation (|6.3|), one obtains from (|8T 



N-l N-l N-l ^_^^| (8.4) 



= / d-ffi-\j)ci-f)^^,E, (^x-E^+J2ij) n (^-7m) + E t{ir,E) n 

The expression in curly brackets is nothing but t(A; E) according to interpolation formula and, 
therefore, the function \E'^ is a solution to the spectral problem ( p.lOl) . Note that the explicit 
form of C(7; E) is not required here. ■ 



^'^ We discard unessential numerical factor in (6.5) 
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Proof of Theorem |3.2| . The proof is heavily based on conjecture that aU the zeros of the Hill 
determinant are simple 0. To derive the formula ( |3.6| ) from ( |3.1j ) let us consider the integral 
over 7i first selecting the expression 

c+(7i) -^c„(7i) 

N - / X (8.5) 



n sinh|(7i-,5, 



k=l 

By construction, it is an entire function but the expressions 

c±(7i) 



N „/ N (8.6) 



n^sinh|(7i-5, 

has an infinite number of poles. Clearly, the original integral can be written as follows while 
selecting the variable 71 and denoting the rest of variables as 7' : 



T 



f J J ' c+(7i) -^c_(7i) _i , 

J d-fidj— fi (7)C(7)^ 

-«3 n sinh J- (^7i - 4) 

'"^^ (8.7) 

/■ r r / C+(7l)/X~^(7) , C f ^ ^ ' g-(7l)/^~^(7) r^f ,^ ,., 

= J d-fid-f — C(7 )^^ + ^ J d-fid-f — C(7 ) ^^ 

c+ n sinh ^(71 - 5k) c_ 11 sinh f (71 - 4) 

k=i "'^ ' k=\ "'^ ' 

where the contour C+ encloses the poles bk + ^^^fc, (^fc >0, A; = 1,...,A^) while the contour 
C_ encloses the poles bk + ^^?T-fc, (^fc <0, A; = 1,...,A^). Such deformations of contours are 
possible since the integrands in the r.h.s of (|8.71 ) vanish exponentially on large upper and lower 
semi-circles respectively; this follows from the asymptotics ( [7.8| ) and ( |7.9| ). Now the integrals 
can be calculated using the residue formula. Note that ^C-{6k + ihnk) = C-^{6k + ihnk) according 
to (|7.13|) . Hence, the final expression can be written in terms of the function c_|_ only. The same 



is true for the rest of integration. 

Let us consider the contribution coming from the residues 5i + ihrii while integrating over 71 . 

Using the explicit expression for the measure 

'jr 

fx-\j) ~ [](7^.-7^)sinh-(7j--7fc) (8 8) 

j>k 

one obtains the following contribution: 

~ Y^ c+{6i+ihni) d'y' JJ {-fk-Si-ihrii) smh.-(-fk - 5i)/i"^(7')C(7')^5i+iani,72,...,7iv-i 

nieZ fc=2 "" 

(8.9) 
where we drop the common factor 



l{N-l)\Y[smh^{6i-6k)\ 



^.10) 



for a moment. Now we consider the analogous integration over variable 72 . It is easy to see 
that there are no poles at 72 = (^i + ihrrik because of the factor sinh |(72— ^1) in the integrand 

^^ Up to now we are unable to prove this conjecture though it seems very plausible. 
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in ( ^.9| ). Quite similarly, while calculating the residues at 72 = ^2 + ihn2 , one concludes 
that the integration over 73 does not include the contribution from the points 61 + ihrrik and 
62 + ihsk and so on. Therefore, the successive evaluation of the residues at "well arranged" 
points 7fe = (5/c + ihuk, {k = 1, . . . , N — 1) gives the contribution 

N-l 

E n C+(4 + ^^^fc)A((5(^) + W^))^5i+inni,...,5^_i+ian^_l (8.11) 

ni,...,njv_ieZ fc=l 

where A(7) is the Vandermonde determinant and the notation y^'^^ means the A^ — 1 di- 
mensional vector obtained from the corresponding N dimensional one by cancellation of the 
s-th component, i.e. y^^^ = {yi, . . . , i/s-i, Vs+i, ■ ■ ■ , Un) ■ After a little algebra one finds that the 
common factor to (^.111) is 



_^^NiN-i)/2 ^ (^_i), -Q sinh^(<5,-4) ) (8.12) 

j>k 



By construction the integrand in (|3.1j ) is symmetric under the permutation of 7 -variables. 
Therefore, one obtains the same answer ( |8.11| ) while calculating all other possible residues from 
the set 8^ ' + ifin^^^ . This gives the additional factor (A^ — 1)! . 

The last step is to take into account the poles at 5jv + ihrif^ . One can easy to show that 
the calculation of the corresponding residues while integrating over 7^ gives the additional 
factor (—1)^"* to compare with ( |8.12| ). Hence, one obtains the final answer (up to unessential 
numerical constant) 



'^j^{xo,x 



E' ^' ' nsinhf(5,-<5,) 
j>k 



N 

xEi-^f'' E A(5(^)+z/inW)C+(5(^)+zW^))xl/,(.)^,,„(.)^^^(xo,a;) 
s=l n(')ez^-i 



(8.13) 



where 



N-l 

C+h) = n c+(7fc) (8.14) 



k=l 



and Theorem 3.2 is proved 
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